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Abstract 

The recent investigations of pure Landau gauge SU(3) Yang-Mills theories which are based on 
the truncated Schwinger-Dyson equations (SDE) indicate an infrared power law behavior of the 
gluon and the ghost propagators. It has been shown that the gluon propagator vanishes (or finite) 
in the infrared limit, while the ghost propagator is more singular than a massless pole, and also that 
there exists an infrared fixed point of the running gauge coupling. In this paper we reexamine this 
picture by means of the exact (non-perturbative) renormalization group (ERG) equations under 
some approximation scheme, in which we treat not only two point functions but also four point 
vertices in the effective average action with retaining their momentum dependence. Then it is 
shown that the gluon and the ghost propagators with the infrared power law behavior are obtained 
as an attractive solution starting from rather arbitrary ultraviolet bare actions. Here it is found to 
be crucial to include the momentum dependent four point vertices in the ERG framework, since 
otherwise the RG flows diverge at finite scales. The features of the ERG analyses in comparison 
with the SDE are also discussed. 
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I. INTRODUCTION 



In quantum field theories, the Green functions are the most fundamental objects in order 
to understand the dynamics. Especially in QCD, the two point functions tell us various 
informations about non-perturbative phenomena, such as dynamical mass generation of 
chiral fermions and also color confinement, which are expected to occur in the infrared 
region. There have been many studies concerning to the dynamical chiral symmetry breaking 
(D%SB) in QCD and its effective models using SDE and ERG. Though these are successful in 
describing D^SB, much less has been well known as for color confinement due to complicated 
gluon self interactions . 

The recent SDE studies for the gluon and the ghost propagators in Landau gauge Yang- 
Mills theoriesjl], 0, 0, 0,0,0,0,0,01 initiated by Smekal, Hauck and Alkofer[l|, indicate the 
power law behavior: the gluon propagator vanishes (or finite) in the infrared limit, while 
the ghost propagator is more singular than a massless pole. IR asymptotic forms of the 
propagators are found to be described by one exponent (ghost anomalous dimension) k as 

Dz{p 2 )-{p 2 Y l+2 \ D G (j?) ~ (p*)- 1 -" , (1) 

for the gluon and the ghost propagators respectively. The actual value of k obtained by 
solving the SDE is dependent on the approximation scheme. However fit of the data by lattice 
simulations at some finite Euclidean momentum region [6] as well as axiomatic consideration 
0] indicates k ~ 0.5. It is also interesting that the solutions of the SDE show consistent 
behavior with the lattice data. Implication of this behavior to confinement is attributed to 
realization of the so-called Kugo-Ojima confinement scenario via the BRS quartet mechanism 
of colored particles [lOj, . In the Landau gauge, the confinement criterion is understood 
as infrared enhancement of the ghost propagator, and has been proven to be fulfilled under 
general ansatz by using the SDE for ghost in the Landau gauge In addition to the 



confinement criterion, the IR asymptotic power solutions satisfy also the Gribov-Zwanziger 
horizon condition [l3l The horizon condition is a dynamical consequence obtained in 
restricting the support of the Faddeev-Popov measure to the interior of the Gribov region 
to avoid gauge copies. The conditions are represented in terms of the propagators as 

lim D z (p 2 ) = 0, lim Ip'DGip 2 )}' 1 = 0, (2) 

which are found to be fulfilled by the SDE solutions. Thus, one may say that the power 
solutions of the SDE most probably capture an aspect of confinement [l5| , though the string 
tension has not been derived directly from the two point functions. 

Here it would be worthwhile to mention also Zwanziger's observations that the cutoff 
prescription at the Gribov horizon resolves an ambiguity in the solutions of SDE 0, 0]. 
First it is noted that cutting off the functional integral at the Gribov horizon does not alter 
the SDE. Also the boundary contribution does not affect the form of SDE, because the 
Faddeev-Popov measure vanishes on the boundary. However, it is shown that the cutoff 
at the Gribov horizon provides supplementary conditions for the solutions of SDE so as to 
make the Faddeev-Popov measure positive definite. Indeed it is needed to know the infrared 
asymptotic power law behavior of the solutions in the practical analyses of the SDE. However 
there are found multiple solutions satisfying the power law ansatz. Then Zwanziger pointed 
out that infrared asymptotic solution in the SDE analyses done so far has been implicitly 
chosen to satisfy the supplementary condition (3,0,111. Thus, in the practical SDE analyses, it 
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is important to impose proper infrared boundary conditions in order to find out the physical 
solution j^jj. Contrary to this, however, the exact renormalization group (ERG) approach 
does not need to use the supplementary conditions for the infrared behaviors, the power law 
solutions with a certain exponent emerge uniquely as will be shown in this paper. 

The ERG describes continuous evolution of the coarse grained effective action in the 
spirit of Wilson's renormalization group^6|. In continuum field theories, this evolution is 
most conveniently expressed in terms of the effective average action T\ with the coarse 



grained (infrared cutoff) scale A [17|,|18(. in practice the ERG has been frequently used in 
various non-perturbative calculations and found to be rather useful. In contrast to the SDE 
framework, in the ERG, solving the flow equation starting with a ultraviolet bare action 
Ta gives a unique infrared solution Fa=o as Ao —>■ 0. Indeed it can be shown schematically 
that the solutions of ERG satisfy the corresponding SDE on the exact generating functional 



level |20Ll21|. However the correspondence becomes totally unclear, once some approximation 



or truncation is performed, except for some rare cases 23|. Especially in the cases that the 
SDE allows multiple solutions, it would be an interesting problem to see whether the physical 
solution can be obtained by solving the ERG equations. Another good feature of the ERG 
scheme is the infrared attractive property of the RG flows. In renormalizable quantum 
field theories, the same infrared effective action is derived from rather arbitrary ultraviolet 
bare actions owing to the universality argument j23[. Thus it is expected to be rather 
straightforward to analyze the infrared two point functions of the Yang-Mills theories in the 
ERG framework. 

In this paper, we investigate the infrared behavior of the effective action, specially the 
gluon and the ghost propagators of Landau gauge Yang-Mills theory by solving the approx- 
imated ERG equation. The main purpose of this paper is to see, first whether the infrared 
power behavior can emerge also in the ERG framework at all, and then which kind of cor- 
rections are necessary to be taken into account for that. Previously, Ellwanger, Hirsch and 
Weber also have studied infrared behavior of the propagators by using ERG |20|, and they 
concluded that the confining 1/p 4 behavior for the heavy quark potential which was derived 
from the effective four quark interactions. However, their RG flows become singular at a 
finite scale, and therefore the infrared power law behavior has not been seen. The most 
different point with our analysis is inclusion of the four point vertices among gluons and 
ghosts, which are necessary ingredient to produce the power law solutions. In their analysis, 
only the four gluon vertex, which was determined by using the Slavnov- Taylor identities, 
was included. Our analysis indicates that the effective four point vertices generated by 
box diagrams play a specially important role, because the large ghost anomalous dimen- 
sion enhances these four point vertices. Indeed, the infrared solutions with power behavior 
are found to be obtained by taking these into account, although with a small exponent 
k « 0.146. 

Application of the ERG to Yang-Mills theories immediately faces problems of the gauge 
or the BRS invariance, because the infrared cutoff of momentum breaks the local symmetry. 
Schematically this problem may be managed by use of the modified Slavnov- Taylor identities 
(mSTi) [H, HQ, which guarantee the broken BRS symmetry of the effective average 



action Ta at an intermediate infrared cutoff scale A ^ to recover the symmetry at infrared 
limit A — > 0. There the BRS non- invariant counter terms, such as gluon mass, should be 
controlled by the mSTi. However, the infrared BRS invariance is not maintained any more 
once we truncate the effective action. Effectiveness of the mSTi in the practical sense is 
unclear apart from perturbative analyses. Therefore, in this paper we do not pursuit for the 
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problems of gauge invariance and simply discard the BRS non-invariant corrections. This 
is a difficult aspect of non-perturbative analysis of gauge theories not only in the ERG but 
also in the SDE framework. 

Furthermore, for the sake of simplicity, we neglect also corrections with the three gluon 
vertex in the flow equations, though the resultant flows will no longer be Yang-Mills theory's 
ones. However, according to the previous SDE analyses 0, 0, IE HI 0|, contribution of the 
three gluon vertex was found to be sub-leading as far as the IR asymptotic solutions are 
concerned. Therefore we also adopt such a simple approximation as the first trial to apply 
the ERG equation for the infrared dynamics of Yang-Mills theories. 

The paper is organized as follows: In section 2, we explain briefly the exact flow equation 
for Landau gauge Yang-Mills theory and also introduce some basic notations. In section 3, 
we present our approximation scheme and give the practical RG equations, which we shall 
analyze in the followings in detail. Results of the numerical analysis are presented in section 
4. There the RG flows of the momentum dependent gluon and ghost form factors and also 
the running gauge coupling obtained in our scheme are shown. Section 5. is devoted for 
some discussions and remarks on further issues. Some comparisons between the ERG and 
the SDE formalisms are also summarized there. 



II. EXACT RENORMALIZATION GROUP 

The exact renormalization group is a realization of the Wilson renormalization group 
transformation in continuum field theories, and there have been known various formulations. 
The most commonly used form is the so-called flow equation written in terms of the effective 
average action TaIhI]. The effective average action T\ is the effective action T obtained after 
only quantum fluctuations with momentum scale p 2 > A 2 are integrated out. Therefore, it 
is an efficient tool interpolating the classical bare action S and the quantum effective action 

r; 

S bare [$] , limr[*;A]=r[*]. (3) 



lim r[$;Al 



Here dependence on the ultraviolet cutoff A is suppressed implicitly. Explicitly the effective 
average action T\ is defined through the Legendre transformation of the IR regularized 
generating functional for the connected Green functions W\; 



W[J; A] = log J V® exp {-£[*] - AS[$; A] + J ■ <&} 



(4) 



where AS^*!?; A] is the IR cutoff term introduced to suppress the infrared mode of p 2 < A 2 , 
and which is normally quadratic in fields The Legendre transformation of W\ is defined 
by 

r[*;A] = -W[J;A] + J • *- A5[*;A]. (5) 

The flow equation is given by infinitesimal variation of T\ with respect to the IR cutoff scale 
A and it can be written down as one-loop exact form. In the case of pure Yang- Mills theory, 
the flow equation may be given as 



9 A r[$;A] = istr 



dR 



5 2 T 



A] + R 



(6) 
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where we introduced the matrix notation following Ellwanger et al. 24]|. The supertrace 
includes momentum integration as well as summation of color and Lorentz indices. The 
field <fr and <fr denote the gluon, ghost and anti-ghost fields in short hand notation as 



$ = {A, c, c) , $ = (A, -c, c) 



(7) 



The symbol R stands for the diagonal matrix of the cutoff functions, which are explicitly 
given by 

R(p) = diag (i#(p; A) , -i^fa A) , -i?*(p; A)) , (8) 
0R(p) = diag (d A Rf v (p; A) , -d A iT 6 (p; A) , -d A R ab (p; A)) . (9) 

The cutoff term added in the definition of the cutoff generating functional given by eq. (jlj 
is written down in terms of these functions as 



AS[$;A] 



1 



Al{-p)R%{p-k)A» v {p) - I c a (-p)R ab (p;A)c b (p) 



(10) 



Here we adopt the cutoff functions of the following form; 



*£(p;A) 



PyPu 



V 

2\J2 cab 



Mp 2 ) + ms RA (p 



p2 



2 xab 



2. r 
p 



R ab (p;A) = R A (p 2 )p 2 6' 



(11) 
(12) 



for the gluon and the ghost fields respectively. We may choose a general cutoff functions for 
i?A here, however in the following analyses, we will take sharp cutoff limit given by, 



Ra(p 2 



lim 

a— >oo 



P 
A2 



(13) 



for the calculational simplicity. 



III. APPROXIMATION SCHEME 

In the practical analysis, we need to perform some approximation in order to solve the 
flow equation schematically given by ©. In the ERG framework it is done typically by 
truncating interactions in I\ at some finite order, which should be compared with truncation 
of equations in the SDE formalism. We may improve the approximation systematically by 
increasing the interactions to be taken in. The recent analyses of the SDE treat only the 
two point functions of the gluon and ghost fields as dynamical ones. However it does not 
mean that we may truncate the effective action or corrections to the action at the two point 
function order. Although the solutions of the ERG and the SDE should coincide to each 
other at the full generating functional level, the truncation ruins the correspondence. As 
we will show in the subsequent section, it is necessary to include the four point vertices to 
realize the infrared finite power solution in the ERG formalism. 

One of our main interests is to see how the infrared power behavior of the gluon and ghost 
propagators emerge in the ERG framework. We shall now explain briefly the approximation 
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scheme which we adopted for this purpose. Our approximation scheme is, in short, to treat 
all diagrams that contribute to the gluon and ghost two point functions. 

Because the flow equation © is one loop exact, the vertices contributing to the corrections 
of the two point functions are not more than four point. Therefore, we truncate I\ at the 
forth power of fields and the explicit form is given as follows, 

r[A,c,c;A] = ~jf r«t(p;A)A«(-p)4(p) 

+ / (27r)V(p + q + r) V$ c p (p, ?. r '> A ) A^p) A b v {q) A c p {r) 

Jp,q,r 

+ f (2vr) 4 5 4 (p + q + r + s) V$f a (p, q, r, *; A) A«(p)A b u (q)A c p (r)A d a ( S ) 

Jp,q,r,s 

- [ f ab (p;A)?(-p)c b (p) 
Jp 

- I (2n)*5(p + q + r) T^(p\q, r; A) A»c b (g)c c (r) 

Jp,q,r 

+ f (2vr) 4 5 4 (p + q + r + s) T« bcd (p, q\r, s; A) A«(p)A b u (q)c c (r)c d (s) 

Jp,q,r,s 

+ J (2tt) 4 5 4 (p + q + r + s) T abcd {p, q, r, s; A) c a {p)c b {q)c c {r)c d {s) 

«™ (t^/i&-7(&). (") 

where the vertices are defined (anti)symmetric, e.g V® b <£(j>, q, r, s; A) = V b ^ d (q,p, r, s; A) = 
• • •, and so on. Substituting the action ((H)) to the flow equation (JHJ) and expanding the RHS 
of the flow equation in powers of fields, then we obtain the coupled flow equations for the 
momentum dependent vertex functions. 

However we may simplify the equations by using more specific form of the two point 
functions. For the gluon two point function, we choose, 



rife A) 



p 2 5 ab . (15) 



Although we write also the A-dependent gauge parameter £(A), we will consider only the 
Landau gauge limit case £(A) — > in this paper. It is noted that this limit can be taken 
consistently, because the Landau gauge is a fixed point f*(A) = 00 under the RG 
evolution. 

To be more precise, the equation (|T3|) is incomplete even in the Landau gauge, because the 
cutoff term (Jllj) induces the gluon mass term ~ m(A) 2 A a A a in the averaged effective action 
due to lack of manifest BRS invariance. Such BRS non-invariant terms must be controlled by 
the mSTi consistently with the RG evolution [13]. Here, however, we neglect simply such BRS 
non-invariant terms and also the corrections to them. Then by contracting the flow equation 
for the gluon two point function with Brown-Pennington tensor B^ v {p) = (5^ — ^PuPv/p 2 ), 
we may project out correction to the gluon mass term. For the ghost two point function, 
we choose, 

f ab (p;A) = 5 ab p 2 f G (p 2 ;A). (16) 

Thus the two point functions are reduced to A dependent one component functions fz and 
fa. Integrating the flow equations toward A — > 0, the gluon and the ghost propagators, 



6 



which are of our present interest, are found to be 

„2\ ^fivip) n ,2\ 



DzMP 



pVz(p 2 ;A^0) 



D G (P 2 



(17) 



where T^ u {p) = (8^ — PfjPv/p 2 ) is the transverse projection. 

Next we restrict the three gluon vertex and the three ghost-gluon vertex to the following 
forms; 

V$ c p (p,q,r;A) = 0, (18) 
T« bc (p\q,r;A) = ig(A)f ab %, (19) 

where the three gluon vertex is just neglected. It is because the gluon loop corrections give 
sub-leading contributions concerning to the infrared momentum region where the power 
behavior is seen in the SDE 0, 0, H 0, HI • Of course the three gluon vertex is not negligible 
at ultraviolet and intermediate momentum region, and therefore resultant RG flows cannot 
be regarded as those for the Yang-Mills theories. However we examine the ERG by imposing 
this condition just for simplicity as the first step of the analyses. While the three ghost- 
gluon vertex is given by the bare form ([190. Owing to the non-renormalization theorem 
for the three ghost-gluon vertex in Landau gauge [25^. the coupling constant is not modified 
g(A) = g(Ao), In fact, if we truncate the bare form ([190. the coupling non-renormalization 
d A g(A) = is also shown by using the flow equation similarly to the perturbative one- loop 
argument. 

Now let us derive the flow equations for the two point functions. From now on we will 
consider SU(3) Yang-Mills theory, although the extension to general SU(N) gauge group can 
be done without difficulty Contracting the flow equation for the gluon two point function 
with BP tensor B^lp), we can easily obtain the following flow equation, 

U ab p 2 d A f z {p 2 -A) 



-^V(Ao) J k k 2 d A R A (k 2 )P G (k 2 ;A)P G ((k + p) 2 ;A) [k a (k + p)pB a p[p)] 
l -5 ab g 2 {A Q ) J k (k + p) 2 d A R A ((k + p) 2 ) P 2 ((k + p) 2 ;A)P G (k 2 ;A) [k a (k + p) P B aP {p)] 
l -J k k 2 d A R A {k 2 )P 2 G (k 2 -A) 
l -J k k 2 d A R A {k 2 )P 2 z (k 2 -A) 



B aP (p)T$ cd (p,-p\k,-k;A)5< _ 



x 12, 



(20) 



for the ghost two point function. The flow equation for the gluon two point function is also 
found to be 

-5 ab p 2 d A f G (p 2 ;A) 

= -36 ab g 2 (A ) ^ k 2 d A R A (k 2 )P 2 (k 2 ;A)P z ((k + p) 2 ;A) [k^T^k + p)\ 

-35 ab g 2 (A ) J (k + p) 2 d A R A ((k + pf) P 2 ((k + p) 2 ; A) P g (k 2 - A) [k,k u T, v (k + p)\ 

«'<•* ( P) -p f k , -k; A) 5^ 



k 2 d A R A (k 2 )P 2 (k 2 ;A) 
\j k k 2 d A R A {k 2 )P 2 z (k 2 -A) 



x 4 



T^(kmr(k,-k\p,-p;A)5 cd x2, 



(21) 
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k 

(-p,t,/3) (p,s,o) (-p,t,/3) //" ^\ (p,s,a) 

WVAAAA < V\/\#( ^/W + ( k 





k+p 



(-p,t) (p,s) (_P,t l-i( ft i P ' S) 





k+p " * 

FIG. 1: The diagrammatic representation of the beta functions for the gluon and the ghost two 
point functions; fz(p 2 ] A) (upper diagram) and fc(p', A) (lower diagram). The crossed circle implies 
insertion of the cutoff function 8\Ra- 

where we introduced short hand notations for the cutoff propagators; 

P '^ ;A ) = WOT (22) 
Pc ( p - A ) ^W;A) + W (23) 

The integrals in the beta functions may be interpreted as specific one loop corrections and 
their diagrammatic representations are presented in FigQ 

Even though the three point vertices are simple in the Landau gauge, it is still difficult 
to solve the flow equations. Because the two point functions and the four point functions 
appearing in the beta functions are fully momentum dependent. In oder to integrate the flow 
equations toward infrared numerically, it is quite preferable to perform the shell momentum 
integrals in the RHS of the flow equations analytically. Then the numerical analyses are 
much simplified. For this purpose, we adopt one dimensional (y-max) approximation and 
take sharp cutoff limit as follows. 

The one dimensional approximation is often used also in SDE analyses, which replaces 
the each form factor as, 

/ 4 ((fc + p) 2 ;A) ^/.(max^V^A), (i = G,Z), (24) 

in the above flow equations (|20|) and (|21|). Next we consider to perform the angular integrals 
of internal momentum. It is noted that the cutoff propagators Pi ((k + p) 2 ; A) are reduced 
to the factorized form as 

P ' (<* + ^ A ) - W /,(ml{^ };A / ' ( {k + P)2 - ^) ' (25) 

in the sharp cutoff limit given by (|13j). The superscript e indicates the sharpness of the cutoff 
function (e ~ a -1 for (|13|)). However, naive sharp cutoff limit induces non-analyticity at the 
origin of momentum space, reflecting the non-local behavior in position space. One way out 
this problem is to use momentum scale p ~ ^/P^Pn instead of the momentum components 
P/i|27(. Here, however, we avoid this problem by simply replacing the all momenta associated 



S 



with the cutoff propagators k + p± + • • • + p n (k is inner-loop momentum, and pi is exter- 
nal momentum) with the inner-loop momentum k. Resultantly the above manipulations 
correspond to replacement of the cutoff propagators as 

^■((^^A)- (<: + riV , (m L {F , p , };A) ^y-A 2 )- (») 

The A derivative of the cutoff propagators are also replaced by 

(k + p) 2 d A R A ((k + p) 2 ) P % ({k + p) 2 ; A) 

1 1 



A'=A 



(k + py[f l ((k + p)i;A) + R A ((k + P y)] 
p)^(max{F,p2} ;A )^ (*' ~ A ') • W 



After this replacement, the shell integrals may be carried out analytically by using the 
formula |27T |: 

5 e (k 2 -A 2 ) h(6 e (k 2 -A 2 )) ^5(k 2 -A 2 ) £ dth(t) as e — 0, (28) 

where h(8 e (k 2 — A 2 )) is an arbitrary function of the step function 9 e . Indeed it is far from 
obvious whether the one dimensional approximation and the sharp cutoff limit introduced 
in our approximation scheme is effective well or not. However it would be inevitable to 
use of these approximations in order to solve the flow equations with retaining the full 
momentum dependence of the two and four point vertex functions. Needless to say this 
approximation may be too brute, although the analysis of the flow equations becomes easy. 
Another approach base on derivative expansion of the effective average action may be also 
possible, however, we would like to leave it for future work pa]. 

Lastly let us consider the four point vertices. By substituting the truncated action (f*HJ) 
into the flow equation (|SjL we may obtain the beta functions for the four point vertices: 
T^ cd (p,q\r,s; A), V^f a (p, q,r, s; A),T abcd (p, q,r, s; A). However it is still difficult to solve 
these flow equations with retaining all the informations of the three independent momenta, 
color and Lorentz indices. Therefore we shall further restrict their functional forms below. As 
stated at the beginning of this section, our truncation scheme deals with only the diagrams 
which contribute to the two point functions. Looking at the flow equations (J2~U|) and (|21|). 
we observe that the four types of four point vertices appear in the flow equations with two 
independent momenta (the external momentum p and the internal momentum k) associated 
with the following momentum channels; 

T$ cd (p, - P \k, -k; A), Vjg(p, -p, k, —k; A), (29) 
T abcd (p, -p, k, -k- A), T$*{k, -k\p, -p; A). (30) 

Here, we have distinguished the internal momentum k from the external momentum p, 
because these momenta are treated asymmetrically in the followings. Then we shall ap- 
proximate the beta functions of the four point vertices by restricting to the above types of 
vertices; we neglect differences in the other momentum channels. Also we restrict the color 
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and Lorentz index structures with the following forms; 





-p\k, 


—k; A) = 


'5 ab 5 cd /8 


V/jtvaP {Pi 


-P,k, 


—k; A) = 


'6 ab 5 cd /8 


T abcd (p, 


-p,k, 


-*;A) = 


'5 ab 5 cd /8 




-k\p, 


-p;A) = 


'5 ab 5 cd /8 



(31) 

(p 2 ,fc 2 ;A), (32) 
(33) 
(34) 

where T^ u (p) denotes the transverse projector and the normalization factors are chosen for 
the later conveniences. Here we introduced the four types of four point vertices Wi, W2, 
W3, W4 as functions dependent on two momentum variables p 2 and k 2 as well as the cutoff 
scale A. Thus we have reduced the four point vertices to these scalar functions by restricting 
structures of the four point vertices, and also by using the one dimensional approximation 
in the sharp cutoff limit. Then the flow equations for the four point vertices are reduced 
to a set of coupled partial differential equations with respect to the momentum dependent 
functions: fc, fz, W±, W 2 , W 3 , W4. One may wonder if the functions W\ and W4 are not 
independent mutually, as is seen from their definition given by (J31|) and (J34)) . However we 
will treat them as independent functions because of the following reason in approximation. 
In order to write down the beta functions for Wi, it is necessary to evaluate the box diagrams 
which consist of the three ghost-gluon vertices (see FigJ2J). Then one encounters the two 
angle integrations which are induced by the propagtors associated with the momenta 
(k' + p) 2 and (k' + k) 2 , where k' is the internal loop momentum in the box diagrams. This 
angular integrals in the box diagrams cannot be performed analytically. Here, therefore, we 
will approximate the integrand by expanding the inner products of momenta that is, 
(k'k) and (pk) to the second order, and then average over the direction of k^. Thus W\ and 
W4 receive different radiative corrections in this evaluation of beta functions, so we treat 
them as independent functions. 

Resultantly, our ERG equations are reduced to a set of coupled partial differential equa- 
tions for the two propagator functions fz(p 2 ;A), fc(p 2 ;A) and four types of the scalar 
functions Wi(p 2 , k 2 ; A), (i = 1, • • • , 4). Diagrammatic representations of the beta functions 
for the two and four point functions are depicted in FigQ and Fig|2] respectively. Concrete 
expressions of the flow equations are summarized in Appendix.A. We will solve this set of 
ERG equations numerically in the next section. 

Before closing this section it is also helpful to define more primitive approximation scheme 
for later discussions. We define the scheme in which the four point vertices are neglected by 
putting Wi(p 2 , k 2 \ A) = under the RG evolution. This scheme is equivalent to take into 
account only the first diagrams of each flow equations for fz(p 2 ]A) and fc(p 2 ',A) given in 
Fig^ At first sight this scheme is more similar to the SDE approximation scheme, which 
has the infrared power solution. However, as we will show in the next section, this scheme 
will fail to reach deep infrared region. The reason will be considered in the last section. 



IV. RESULTS OF NUMERICAL ANALYSIS 



In this section we will present the results obtained by numerically integrating the flow 
equations, specially with paying much attention on infrared behavior of the gluon and the 
ghost propagator functions, fz(p 2 ]A) and /g(p 2 ;A). First of all, we need to specify the 
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P(k) 

Wi (WO 





v * 






® i + 18 i - $> t 





FIG. 2: Diagrammatic representations of the beta functions for the four point vertex functions: the 
4-ghost-gluon vertices W±(p 2 , k 2 ; A) and Wi(p 2 , k 2 ;A) (upper line), the 4-gluon vertex Wi{p 2 ■, k 2 ; A) 
(center line) and the 4-ghost vertex W3 (p 2 ,k 2 ;A) (lower line) . Here p denotes the external and k 
denotes the internal momenta when they are used in the flow equations for the two point functions 
shown in Fig^ 



initial condition for the effective action in order to proceed the numerical integration. If one 
choose a very large scale as the initial cutoff scale Ao, then the boundary condition for the 
flow equations may well be given by the bare action. This is because the arbitrary irrelevant 
terms do not affect the infrared behavior according to the discussion given in section 
Obviously the large initial cutoff scale A is characterized by the small coupling constant in 
asymptotically free theories. We thus define the initial conditions at the initial cutoff scale 
A simply by 

<7(Ao) = 2.0, (35) 
f z (p 2 ;A ) = f G (p 2 ;A ) = l, (36) 
Wi(p 2 ,k 2 ;A ) = (i = l,-..,4). (37) 

Here it should be noted also that the bare four gluon vertex does not appear in the ini- 
tial conditions (J37|) . The presence of such a momentum independent four point vertex is 
irrelevant in the flow equations, because the corrections are contracted by the B-P tensor 

bm. 

It would be better to use the initial conditions defined with including quantum fluctuation 
from larger momentum region, p 2 > Aq. This is easily carried out by solving the perturbative 
flow equations under suitable renormalization conditions [20L l26| . However such corrections 
are relevant for flows only at the ultraviolet region. Our main concern is to see whether the 
infrared power solution is realized as an infrared attractive solution in the ERG framework. 
Therefore we employ the above simple initial conditions (|35)) M|37 )) . 

Now let us present the results of numerical analysis of the flow equations. In Figs|3]and|U 
we show flows of the gluon function fz(p 2 ', A) and the ghost function fcip 2 ] A) respectively 
with varying the infrared cutoff scale A = exp(— £)Aq. In the figures their snapshots at 
t = 0.00,3.12,4.16,4.94,8.06,13.00 are presented in one panel. From these flows one may 
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observe that the propagator functions fz(p 2 \A) and fG{p 2 ',A.) freeze immediately once the 
cutoff scale is lowered than the momentum scale, namely, at the region of A 2 < p 2 < Aq. 
The graphs are strongly bended around p 2 ~ A 2 and become almost flat at lower momentum 
region. It may be supposed that the bends are due to the crude approximation. Because the 
integrand of these beta functions are infrared cutoff, and therefore the corrections around 
the p 2 ~ A 2 may be overestimated due to the one dimensional approximation in the sharp 
cutoff limit employed in the previous section. Anyway, as lowering the cutoff scale A, it is 
found that the frozen solutions drastically change around p 2 /Aq ~ 0.0001 and eventually 
show the power scaling. This means also that power behavior is certainly realized in the 
infrared propagators defined by (fTTj) . 

In FigElthe result for these propagator functions obtained at t = 13.00 are shown. There 
we also show these functions evaluated with neglecting the four point vertices Wi(p 2 , k 2 ; A) 
in the flow equations for comparison. It is seen clearly that these functions enjoy power 
behavior in the infrared without singular poles. On the other hand if the four point vertices 
are neglected, we observe that the ghost propagator function diverges at a finite scale, and 
our numerical calculation is broken down. In the present case, the singularity appears at 
t w 5.00. The snapshots in Fig|5] is given adjacent to the singularity. It is seen that 
the two schemes present qualitatively the same behavior at the ultraviolet region, however 
the evolutions in the infrared become drastically different. Thus we may conclude that 
contributions through the four point vertices are significant in the infrared region to achieve 
power behavior with avoiding singularity. 

Interestingly these infrared power solutions are related with the running gauge coupling 
constant, which may be defined naively in the RG framework by 

U 47r/ z (p 2 = 0;A(t))/ 2 (p 2 = 0;A(t))- 1 ] 
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FIG. 4: The flows (A dependence) of the ghost propagator function /c(p 2 ;A) at various cutoff 
scales t = 0.00, 3.12, 4.16, 4.94, 8.06, 13.00. 



The running coupling obtained by the present analysis is shown in Fig El compared with that 
in the scheme neglecting the four vertices. We see that the running coupling converges to 
an infrared fixed point a* ~ 4.70. In contrast to this the running coupling without four 
vertex contributions is found to hit the Landau pole. In the SDE approach, the running 
gauge coupling is sometimes defined as a function of momentum, which is regarded as the 
renormalization scale 0,0). In our ERG framework, the coupling constant would correspond 
to 

/ 2\ 9 (Ap) /„„x 

a{P ] ~ ^fz{v 2 ; A - 0)f G (p2; A - 0) ' {6J> 

It is seen also that this coupling constant shows the qualitatively same behavior with a(t) 
and that an infrared fixed point exists. 

From the definition of the running coupling (|3*H|) , it is clear that appearance of the infrared 
fixed point can be realized by exact cancellation of the corrections to the gluon and ghost 
propagators. The infrared power behavior implies that the exponents must satisfy the 
specific ratio as shown in (0). The exponents, that is the anomalous dimensions, may be 
defined also as A-dependent quantity, 



Vz(t) = -~\ogf z (p 2 = 0;A(t)), 



VG(t) = - 1 -^logf G (p 2 = 0;A(t)). (40) 



In the infrared limit (t = 13) these exponents are found to converge to stable constants 
given by 

1 



K = T] G 



-rjz « 0.146, 



(41) 



where k denotes the exponent in Eq. (0) used frequently in the SDE studies . 
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FIG. 5: The propagator functions fz(p 2 ;A) (full line) and /c(p 2 ;A) (long-dashed line) obtained 
at t = 13.00, which show power behavior in the infrared region. For comparison these functions 
evaluated with neglecting the four point vertices are also shown for fz(p 2 ', A) (dashed-dotted line) 
and fc{p 2 ', A) (dashed-dashed-dotted line) just before reaching singularity (t=4.93). 

V. CONCLUSION AND DISCUSSION 

We have performed numerical analysis of the approximated ERG equations for Landau 
gauge SU(3) Yang- Mills theory. There flows of the momentum dependent four point ver- 
tices as well as the momentum dependent gluon and ghost propagators were examined. In 
order to setup the approximated flow equations with retaining momentum dependence of 
vertices in a minimal way, we restricted the diagrams to those which contribute to the two 
point functions. Further we employed one- dimensional approximation and sharp cutoff limit 
to make numerical analysis easier. The resultant flow equations are given in the form of 
coupled partial differential equations which have been solved numerically. We solved the 
flow equations with an ordinary bare action as the initial conditions for the effective action. 
However infrared behavior of the flows are not affected by the arbitrary irrelevant terms in 
the effective action at ultraviolet. 

Our numerical solution exhibits infrared power behavior for the propagators Dz{p 2 ) ~ 
(p 2 )~ 1+2K , Dc{p 2 ) ~ (p 2 ) _1_K with the exponent k m 0.146, which is similar to the result 
obtained by SDE studies. Also the running gauge coupling is found to approach an infrared 
fixed point a* ~ 4.70. Compared with the SDE analyses, the most different point is that 
inclusion of the four point vertices generated during flow is crucial to obtain the power 
solutions with avoiding infrared divergence. In the SDE approach, only the dressed four 
gluon vertex appears through the two loop diagram in the equation for gluon propagator and 
existence of the power solution has been confirmed irrespectively to the vertex contributions 

manias. 

This difference may be considered as follows. The SDEs are a set of integral equations 
among the Green functions, and there all momentum modes are coupled to each other. If 
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FIG. 6: The running coupling a(t) (solid line). The long-dashed line stands for the running coupling 
obtained neglecting four point vertices. 

we solve the SDEs for the propagators iteratively, then the propagotor may be represented 
by an infinite sum of diagrams whose internal momenta are not restricted. On the other 
hand the ERG equation is reduced to a set of the beta functions for the Green functions, 
which are given by one-loop exact forms diagrammatically. Therefore if we could solve the 
flow equations iteratively, then the Green functions are given by an infinite sum of diagrams 
involving all loops, which must be identical to the solution of the SDE. However this is not 
the case after truncation. Note that internal momenta of the one-loop integrals of the beta 
functions are restricted around the cutoff scale A, and contributions from higher momentum 
loops are already taken in the effective vertices in the ERG framework. Also the beta 
functions for the effective vertices involve the higher point vertices as well as the lower ones. 
Therefore, in our present concern, some part of the corrections to the two point functions 
comes through four point effective vertices. If we truncate the effective action and discard 
the four and higher point functions, the corrections generated through higher point vertices 
are lost in the resultant propagators. Speaking more explicitly, the diagrams in which the 
inner loops carry lower momenta than the outer loops are suppressed in the diagrammatic 
terminology. Inclusion of the four point functions improve this situation drastically indeed. 
In the case of Landau gauge Yang-Mills theory, the ghost propagator is considered to be 
enhanced at infrared. Therefore, particularly couplings of the soft ghost modes with other 
high energy modes are expected to give non-negligible contributions. 

Finally we shall discuss residual problems and the advantageous features of the ERG 
approach. Although we found the power solution in the ERG framework, the obtained 
exponent k ~ 0.146 is fairly smaller than the currently most reliable value k ~ 0.5 expected 
from the SDE analysis and also from axiomatic considerations. Furthermore transient from 
the logarithmic to the power behavior of the propagators is not smooth. Probably these are 
due to the fact that our approximation scheme is still poor in the following points. Firstly 
the corrections for the two point functions with p 2 ~ A 2 are fairly influenced by use of the 
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one-dimensional approximation and the sharp cutoff limit. Then the four point vertices 
which have other channels and also higher order vertices may contribute to flows of the 
two point functions. These issues are currently investigated with employing another type of 
approximation scheme jiil]. Further, inclusion of three gluon interactions, which have been 
omitted in this paper just for the simplicity, and also dynamical quarks are remained as 
obvious problems. 

The advantageous features of ERG may be considered as follows. We could obtain the 
infrared power behavior for the gluon and the ghost propagators without any ansatz for their 
infrared behaviors, and by starting from rather arbitrary ultraviolet bare actions. On the 
contrary all the SDE analyses so far were based on the infrared asymptotic analysis, which 
needs an appropriate ansatz. In the case of Landau gauge Yang-Mills theories, the power law 
ansatz was found to suit nicely. However in the other situations such as in the deconfinement 
phase the appropriate ansatz is unknown a priori. It would be another advantage of the 
ERG to obtain some informations of the infrared four point functions at the same time. In 
these respects the ERG formalism seems to be suitable to explore the infrared dynamics in 
more general cases. Especially exploration of dynamics near phase transitions would be quite 
interesting and challenging problem to be considered, and the ERG approach is expected to 
be useful also in such studies. 

Note added : After completion of this paper, Ref.ji^] appeared which addresses related 
issues. 
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APPENDIX A: APPROXIMATED FLOW EQUATIONS 

The flow equations derived by using the approximation scheme explained in sect ion ITTT1 are 
summarized here. The all momentum variables and dimensionful functions are normalized 
by the cutoff scale A, and we use x = p 2 /A 2 ,y = k 2 /A 2 . By contracting the external gluon 
momentum with B-P tensor B^p), all vertex functions Wi(x,y; A) are factorized with the 
momentum variables. So we define the factorized functions by Wt(x, y; A) = xy x Wi(x, y; A) 
and use them below. The flow equations, which are somewhat lengthy though, are found to 
be in the followings. 



d t f z (x;A) = -2xd x f z (x;A) 

o 2 (A ) X( , A w_ v ( 1V jx- 2 (3-2x- 1 ) (x> 1) 



1 ^/^l;A)/G 1 (max(:r,l);A) * j 1 ^ < ^ 

1 ^/G 1 (i;A)w 1 (x,i;A) + ^ a 

d t f G {x;A) = -2xd x f G (x;A) 



f-— fc\l- A) W^x, 1; A) + —J z \^ A) W 2 (x, 1; A), (Al) 



g 2 (A ) , i . . fx" 2 (x > 1) 

"^~2 /g (lj A)/ z (max(a;, 1);A) x < l K < > 
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1 ; fa\l; A) W 3 (x, 1; A) - J-jf z l (l; A) W 4 (x, 1; A), (A2) 



^i(a;,y;A) 

-4V^(x, y; A) - 2xd x W 1 (x, y; A) - 2yd 2/ iy 1 (x, y; A) 

' /J 2 (l; A) W 2 (x, 1; A)W 4 (y, 1; A) + -^/ G 2 (l; A) W^x, 1; A)W 3 (y, 1; A) 



-^/z 2 (l; A)/,'.n,ax ( ,. 1); A) W 2 (x, 1; A) x { ^ ^ > }) 

+^/ G - J (l; Atf^maxfo, 1); A) ^(x, 1; A) x | f * g > J] 

3(? 2 (A ) 2 x - J (3x~ 2 - 2x~ 3 ) (x>l) 

+ 47r2 /g (!; a )/g (max(x, 1); A) W 3 (y, 1; A) x j J ^ < ^ 

+^^/g 2 (1; AJ/o^maxCx, 1); A)/^ 1 (max(y, 1); A) 
/ (l+y(6 + y)) \ f (3x~ 2 - 2x~ 3 ) (x > 1) 

I (i + y) 4 ) U (x<i)' (A3) 

c> t W 2 (x,y;A) 

-4W 2 (x, y; A) - 2x<9^ 2 (x, y; A) - 2y ( 9 J/ jy 2 (x, y; A) 

-^/ z 2 (l; A) iy 2 (x, 1; A)W 2 (y, 1; A) + ^/^(l; A) W^x, 1; A)^(y, 1; A) 

g 2 (A ) 2 .w.i, , ^ \ * \ >7t / , .s ( (3x~ 2 - 2x~ 3 ) (x > 1) 
^t/g (!; A)/ G 1 (max(x, 1); A) W 1 (y, 1; A) x | \ J ^ < ^ 

^/ G - 2 (l;A)/-( m ax(y,l);A)iy l( x,l;A) x { " ^ (j > J) 

,.4 / A \ 



^^/ G 2 (l; A)/ G 1 (max(y, 1); A)/ G 1 (max(x, 1); A) 
/ (5 + 3y)\ / (3x~ 2 - 2x~ 3 ) (x > 1) 

x V(i + y)V x I 1 (^!)' 



d t W 3 (x,y;A) 

-4W 3 (x, y; A) - 2xd x W 3 (x, y; A) - 2yd y W 3 (x, y; A) 



(A4) 



' f z 2 (l; A) iy 4 (x, 1; A)W 4 (y, 1; A) + -^/ G 2 (l; A) W 3 (x, 1; A)iy 3 (y, 1; A) 



48tt 2 ^ v ' ' "*v~>->-/"*v*,-,~y ■ 27r2 
f ^^ 2(1; A )^ 1 ( max ( a; ' ^ A ) ^ 1; A ) X { (x < 1) 



^^/g 2 (1; Vfz\^(y, i); a) w 3 (x, i ; a) x j j g - ^ 

^/z 2 (l; A^Vaxtx, 1); A) iy 4 (y, 1; A) x / " ( * > 1} 



g 47r 2 JZ \-,"JJG \" — VHV "4 W VV - \(3-x) (x<l) 
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" 2(Ao) -/z 2 (l; A)/cVax(*, 1); A) W 4 (y, 1; A) x / ^ ~ ^ <» > 



647r 2 ^ ^V-/G ^x;,^,,^,.,^^^^ 

+^^/g 2 (1; A^maxfo 1); A)/J 1 ( ma x(y l 1); A) 
x / (l + y (6 + y)) \ x / x- 2 (x > 1) 



^ 4 (Ao) 



16tt 



l+y) 4 ; \i (*<i) 

2 /J 2 (l; Aj/c 1 ^, 1); A)f G \mzx(y, 1); A) 



1 ^ / (32T 1 -ar 2 ) (x > 1) , . , 



cW 4 (x,y;A) 

-4iy 4 (x, y; A) - 2xd x W i (x, y; A) - 2yd y W i (x, y; A) 

' f z 2 (l; A) W 4 (x, 1; A)W 2 (y, 1; A) + -^/ G 2 (l; A) W 3 (x, 1; A)^, 1; A) 



•^^/z 2 (l; Al/c'M,, 1); A) l;A)x{ ^ [l < !) 

+^^/ G 2 (l; Aj/^M,, 1); A) ffi(y, 1; A) X | f ' > j } 



32v 

■^ /G - 2 (l;A)/ G -Vax(y,l);A)^l;A) x { for""*-) (y > 1) 
3s 4 (A 



0J / G 2 (l; A)/ G 1 (max(y, 1); A)/^(max(z, 1)? A) 



r'2 



j7T 

(5 + 3y)\ [ x" 2 (x > 1 



HjTTW x \i (x<D- (A6) 
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